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ABSTRACT A positive correlation between the frequency of type-B quasi-periodic oscillations 
(QPOs) and the associated hard power-law flux has been discovered in a few black hole bina- 
ries. The correlation is explained quantitatively based on the Alfven wave oscillation model. 
The QPO is assumed to be generated by the Alfven wave oscillation in the standard thin 
Shakura-Sunyaev disk (SSD) at radius with the maximum radiation flux. The power-law flux 
is calculated based on inverse Compton scattering of soft photons from SSD by a medium of 
hot electrons in the corona or the base of jet. Both the explicit solutions and numerical re- 
sults of the correlation between the QPO frequency and the power-law flux with the changing 
accretion rate are obtained. The modelled correlation fits well into the observed one within 
reasonable parameter range. The positive correlation could be interpreted by the hypothesis 
that a stronger magnetic field leads to a higher Alfven wave frequency and a higher electron 
temperature and therefore a stronger power-law flux. The results suggest that type-B QPOs 
may be related to the toroidal magnetic activities in the accretion disk or the jet. 
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1 Introduction 


Complex spectral and timing properties have 
been observed in Black Hole Binaries (BHBs). 
During an outburst, three stages are classified, 
ie., the hard state, the soft state and the steep 
power-law state!!! or the intermediate statel? 3]. In 
the hard state, the spectrum is dominated by a 
non-thermal power-law component with a typical 
photon index —1.7!!. In the soft state, the spec- 
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trum is commonly described as ~1 keV thermal 
emission and can be well interpreted by the stan- 
dard thin Shakura-Sunyaev Disk (SSD) modell“. 
The spectrum of the intermediate state consists of 
both strong thermal and non-thermal components 
with typical photon index ~2.5 and the luminosity 
is high. The intermediate state is further divided 
into the Hard InterMediate State (HIMS) and the 
Soft InterMediate State (SIMS) by Belloni et al. 
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'The HIMS is spectrally dominated by a hard com- 
ponent and with band-limited noise in the Power 
Density Spectrum (PDS), and the SIMS is spec- 
trally dominated by a soft component and with 
power-law noise in the PDS! 51, These spectral 
states appear in counterclockwise order in the q- 
shaped Hardness-Intensity Diagram (HID) in the 
sequence of “hard state > HIMS 一 SIMS — soft 
state — SIMS > HIMS — hard state” H 3 5-9], 
The hard and 
strong variabilities and quasi-periodic oscillations 


intermediate states show 


(QPOs) are usually observed. A narrow peak in 
the PDS is called a QPO if the quality factor 
Q > 2 (Q is defined as the ratio of the centroid 
frequency and the full width at half-maximum of 
the peak). Low-frequency QPOs (LFQPOs, 0.1- 
30 Hz) are common in the hard and intermedi- 
ate state, while high-frequency QPOs (HFQPOs, 
> 30 Hz) generally appear in the steep power-law 
state!’ 8 and SIMSI. LFQPOs are classified into 
type A, B and C according to the characteristics 
of the PDS and phase lagsl? 11. Type-C QPOs are 
the most common type of LFQPOs with Q 2 8 
and centroid frequencies of 0.01-30 Hz, which are 
observed in the hard state and HIMS"!13!. Type- 
B QPOs are observed in the SIMS with Q > 6 
and typical centroid frequencies of ~6 Hz or 1- 
3 Hzl* 14]. Type-A QPOs are characterized by a 
weak and broad (Q < 3) peak with centroid fre- 
quencies of ~6-8 Hz, which usually appear in the 
soft state and SIMS just after the transition from 
the HIMS!%: 151, 

The physical mechanisms of QPOs are not 
well understood although a variety of models have 
been proposed. It is believed that HFQPOs are 
possibly related to the inner region of the accre- 
tion disk since the period is comparable to the dy- 
namical time scale of the inner disk! 7], Possible 
mechanisms for LFQPOs especially type-C QPOs 


include the relativistic precession"9), the accretion- 


17] global disk oscillations", 


ejection instability! 
the oscillation of shocks in accretion flowsU?!, 
oscillations of the transition layer between the 


P9. the magnetohydrodynamics 


disk and corona 
(MHD) Alfven wave oscillations in the accretion 
disk?4, etc. Ingram et al.??! improved the rela- 
tivistic precession model for type-C QPOs by con- 
sidering the Lense-Thirring precession of a radially 
extended region of the hot inner flow based on the 
truncated disk model, which sucssefully explained 
the frequency range and the spectral behavior of 
type-C QPOsL?22 ?23]. The precession of a small-scale 
jet is also a possible candidate for the production 
of LFQPOs, which can naturally explain the high 
energy, soft lags and energy dependence for fre- 
quency and rms (root mean square) of LFQPOs in 
MAXI J18204-070P4l., 

Type-B and type-A QPOs are poorly under- 
stood. Type-B QPOs are usually observed in the 
transition from the hard to soft state, the appear- 
ance of which indicates the start of the SIMS" 151, 
They usually appear when BHBs reach their max- 
imum luminosities“! and occur almost simultane- 


”[6]. Therefore, the inter- 


ously with the “jet line 
pretation of type-B QPOs are important to the 
study of the states and state transitions of BHBs. 
Type-B QPOs have completely different proper- 
ties with type-C QPOs in many aspects!!+ 25-26], 
and there is growing evidence for the jet origin of 
type-B QPOs?>?7], Motta et al.P5] analysed the 
LFQPO data of 14 BHBs and found that type- 
B QPOs are stronger when the accretion disk is 
closer to face-on, which supports the hypothesis 
that type-B QPOs are related to the jet. Russell 
et al.?9! discovered a compact jet in the SIMS con- 
temporaneous with a type-B QPO in the 2002 out- 
burst of 4U 1543-47. Homan et al.!?9! found that a 


discrete jet ejection launched within ~2—2.5 hours 
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after the appearance of a type-B QPO in MAXI 
J1820+070. Kylafis et al.8°! successfully explained 
the periodic variation of photon index I with the 
type-B QPO frequency in GX 339-4 based on a 
precession jet model. T'he association of type-B 
QPOs and jet is also evident for H 1743-322 and 
MAXI J1348-630P133], 

A strong correlation between the type-B QPO 
frequency and the power-law flux (hereafter “the 
QPO-PL relation”) has been discovered for a few 
BHBs'^ 34, which suggests that type-B QPOs 
may be related to the inverse Compton scat- 
tering either from the base of the jet or the 


ali 25, 34-35]. However, there is a lack of 


coron 
quantitative calculation to interpret this corre- 
lation. In this paper, we attempt to explain 
the QPO-PL relation quantitatively based on the 
MHD Alfven wave oscillation modelP 36739], in 
which the QPO is produced by the toroidal Alfven 
wave oscillation due to the radial perturbation of 
the accreting matter in the inner region of the ac- 
cretion disk, and the “soft” photons from the disk 
are inverse Compton scattered by a medium of hot 
electrons to generate the power-law flux. 

The paper is organized as follows. The de- 
scription of the model is given in Section 2. Then 
the QPO-PL relation is calculated and compared 
to the observed one in Section 3. Finally the con- 


clusion is given in Section 4. 


2 The model 


2.1 The power-law flux 

In the relativistic standard thin disk model, 

the radiation flux per unit disk face area is^?! 
MM 
B ES (1) 
where G is the gravitational constant, M is the 
black hole (BH) mass, M is the disk accretion rate, 
and f is a function of dimensionless disk radius 


Í, 
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r = R/R (Rg = GM/c? is the gravitational ra- 
dius, where c is the speed of light) and BH spin 
a, given in Ref. [40]. The flux F is assumed to be 
radiated as blackbody with temperature 


T(r) = ( ) = ( s) i , (2) 


where ø is the Stefan-Boltzmann constant. The ra- 


3GMM 
Sr Ro 


F 


oO 


diation intensity is 


(3) 


where h is the Planck constant, k is the Boltzmann 


2hv? | 


hw =l 
= 1 kT(r) — 
1, = P feats i], 


constant and 1; is the frequency of the emitted disk 
photons. 

We assume the soft photons emitted from the 
surface of the accretion disk are inverse Compton 
scattered by some hot thermal electrons with tem- 
perature T, in a medium above the disk like corona 
or the base of the jet. If the hot electrons are non- 
relativistic and the soft photons from the disk sat- 
isfy hy, € hvs < kT. (where vs is the maximum 
frequency of the soft photons), the Comptonized 


spectrum has a power-law form!*!*); 


7 -r 
bak (2) 
Vs 


The typical value of the spectral index in the 
SIMS associated with type-B QPOs is LT ~ 1.5 
(corresponding to photon index ~ 2.5)! 5 341, We 


(4) 


take this value in our calculations. 
We take v, as the frequency corresponding to 
the maximum radiation intensity at the hottest 


part of the disk. By Wien's displacement law, we 


have 
3GMMY 
hvs ~ 3kTmax = 3k E 1/4 (a). (5) 
That is 
Vs = 6.55 x 1015m- "^5 ^F1^(a.) Hz, (6) 
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where m = M/Mo (Mg is the solar mass), m = 
M / Mgaa (Mgaa = 1.4 x 105m g.s-! is the Ed- 
dington accretion rate), Tmax is the maximum tem- 
perature of the disk, and finax(a.) is the maximum 
value of f, which is a function of a,. 

'The observed power-law flux is 


n 


FpL 


2hy31T Re cosi(vyl T — vt) 
i c?d?(1—T) 
Tout hus si 
i, Jers = 1 277d7 ， (7) 


where d is the distance to the earth and 7 is the 
inclination of the line of sight to the normal of the 
disk. ri;z Fi, / R4 and ro Rout/ Rz are the inner 
and outer dimensionless radii of the disk respec- 
tively. We take rin as the radius of the innermost 
stable circular orbit rigco and rout = 1000. 


2.2 The QPO frequency 


In the Alfven wave oscillation model, an 
Alfven wave oscillation propagates along the 
toroidal magnetic field lines in the inner accretion 
disk due to the perturbation of radial velocity of 


the accreting matter, and the QPO is produced by 


the Alfven wave oscillation with frequency?! 39] 
Ww k'va 
UE 8 
"REO 20 20 (8) 


where k’ ~ 2x/R is the wave number and R is 

the disk radius at which the oscillation occurs, and 

va = B/\/4rp is Alfven speed, where B is the mag- 

netic field and p is the mass density of the disk. We 

adopt the generalized viscosity prescription? “41: 
2 

B = app , 


where a is the viscous parameter, Pio, = 


(9) 


Peas + Praa 
is the sum of the gas pressure P,,, and the radi- 
ation pressure P444, and 0 < w < 1 is a free pa- 
rameter. The magnetic field B and mass density 


p can be obtained by solving the relativistic ac- 
cretion disk Eqs. (5.6.14) and (5.8.1) in Ref. [45] 
and Eq. (9). Explicit solutions for special cases and 
typical values of u are listed below. 

(1) Radiation pressure dominated region 


(Prot ~ Paa, the Rosseland mean opacity & 


1 


Res = 0.4 cm? - g~!, where Res is the corresponding 


opacity by scattering from free electrons). 
For 1 = 0, we have 


T, = (5.0 x 10" K)o- 4m V/4p-3/8 
_4-U2B1251/4 ， 
p = (1.7 x 1077 g: cm 3)a lim- 
A BDE? L? , 
B = (6.2 x 108 Gs)m7"?r 3/4 A IBE"? , 
Varo = (2.9 x 10° Hz)a!/?m- rr 9/2? 


AB p gie. 


145,-2,3/2 


(10) 


where T, is the central temperature of the disk. A, 

B, D, € and £ are relativistic correction factors 

given in Ref. [45], which are functions of r and a,. 
For u = 0.5, we have 


T, = (8.1 x 108 K)o?/?m ?/9452/9,2/3 
A7 9/9]-1/9g1/9 p2/9 
p = (2.5 x 1074 g.cm am 9g 100 


r1/3 4726/9 g4D5/9¢£13/9 p—10/9 , 


B = (6.2 x 108 Gs)m-1?5—3/^4 4-1 BE! , 
Voro = (7.5 x 104 Hz)o/t/9gy-19/185,5/9,23/12 


A4/9g-1p-5/8g-2/9 75/9 . (11) 
For 1 = 1, we have 


T. = (1.3 x 109 Ka om 71 5552/5,9/19 
B-2/5D-1/5 [2/5 
p = (0.09 g - cm-?)ao*/5g 74/555, 72/5,—3/5 A—2 
BPE Dl/5¢ p -2/5 
B = (6.2 x 108 Gs)m-1?7—3/^ A-1ggV? , 
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vopo = (4039 Hz)o?/5g- 1 /105,1/5,729/20 


B-M5]m-1/10p1/5. (12) 


(2) Gas pressure dominated region (Pio, ~ 
Paas, R = Res = 0.4 em?.g ^L). We have (inde- 
pendent of p) 


T, = (1.8 x 109 K)ao-U5m- 72/59/10 
B-3/5]-M/5 g2/5 
p = (61 g. cm7?)a7 7/105 77/1055,2/5,33/20 
A-183/5D-1/5€1/2 ea/s. 
B = (1.6 x 101? Gs)a1/205, -9/20,2/5,.—51/40 


A-1/251/10p-1/5e1/4 p2/5 
voro = (4039 Hz)a?/ mn 1109 1/5, 99/20 


B-M5p-V10p1/5. (13) 


We assume that the QPO is produced at 
the radius rmax where the radiation flux reaches 
the maximum value. We fit the relationship be- 
tween Tmax and a, based on Eq. (1) and obtain 
Tmax © 1.59rigco 一 0.3a,. For a, = 0.5, substitut- 


ing Tmax into Eqs. (10)-(12), we have 

voro = (2373 Hz)a!"?m^!m, (for y= 0), 

voro = (536 Hz)o*/9m-19/155,5/9 | (for p = 0.5), 

voro = (163 Hz)o?/5m-H/105,1/5 , (for u= 1). 
(14) 


The numerical results of B, p and vepo at 
Tmax as functions of m with different parameters 
(m, a,, œ and yw) are shown in Fig. 1. We see 
that B and vepo are positively correlated with 
m, while p changes nonmonotonically with rn, first 
increases with m, and then decreases after reach- 
ing the maximum where the radiation pressure be- 
gins to dominate over the gas pressure. B, p and 
YVepo are all anti-correlated with m and positively 
correlated with a,. B is insensitive to a and p. p 
is anti-correlated with a and positively correlated 
with u. vapo is positively correlated with o and 
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anti-correlated with u. The slope of the vopo-m 
relation is steeper with a smaller p. 


3 The QPO-PL relation 

In order to compare our results to the ob- 
served QPO-PL relation in Ref. [34], we adopt the 
same energy range (3-200 keV) for the power-law 
flux. Substituting vı = 3 keV/h, və = 200 keV /h, 
T = 1.5 and vs in Eq. (6) into Eq. (7), we obtain 


“a ea.) 


1000. = 
f Jers — : 27r7d7 . 


TISCO 
Combining Eqs. (14)-(15), we have the Fpr- 
vopo relation for different values of r and y as fol- 


For, = 3.13 x 1099 7/5,5/5 208? 


(15) 


lows (assuming a. = 0.5, i = 50? and d = 5 kpc). 
(1) m Z 0.1, p=0: 


2 y?/8 Ot V 


-9/6,, 
/ Qeo TE: Cm ^. s 


Fp, = 6.5 x 107a 


(2) m > 0.1, u = 0.5: 


Fe, —12x 107 139,79/10,,,241/80 


81/40 E 


=i 
Vopo erg:cm “Ss . 


(17) 
(3) m2 0.1, p=1 or 107 $m S 0.1: 


Fo, =1.5x 10729 qy-9/4 777 113/16 


45/8 mm 


Vopo erg .cm ^. S (18) 


We can see that the slope of the QPO-PL rela- 
tion becomes steeper with the decreasing accretion 
rate. Our model predicts Fp, c Daas with mod- 
erate accretion rate. 

Numerical results of the modelled QPO-PL re- 
lation for different values of m, a,, Q and p are 
shown in Fig. 2. We find that the influence of a, 
and o on the relation is similar. Both larger val- 
ues of a, and a cause the curve move to the right 
of the graph, while a larger m value causes the 
curve move up and left. A larger u value leads to a 
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steeper correlation when the accretion rate is high. dependent of u when the pressure is dominated by 
Three curves with different values of u coincide the gas pressure as shown by Eqs. (13)-(14). 


when the accretion rate is low because vapo is in- 


08 
.00 0.02 0.04 0.06 0.08 0.10 0.00 0.02 0.04 0.06 0.08 0.10 


m m 


(à) 


B/(10°Gs) 


pK(g:cm 3) 


oOooooooo: 


CON A Oo OOSNORBROO-^ 


Yopo/Hz 


4 


.00 0.02 0.04 0.06 0.08 0.10 


.00 0.02 0.04 0.06 0.08 0.10 

m m 
Fig. 1 The curves of B, p and vapo as functions of m. Panel (a): the dashed, solid and dotted lines correspond to m = 5, 10, 
and 15 respectively with fixed a, = 0.5, a = 0.15 and u = 0.5; Panel (b): the long dashed, dot-dashed, dashed, solid and dotted 
lines correspond to a, = —0.9, —0.5, 0, 0.5, and 0.9 respectively with fixed m = 10, a = 0.15 and u = 0.5; Panel (c): the dashed, 


solid and dotted lines correspond to a = 0.01, 0.15 and 0.3 respectively with fixed m = 10, a. = 0.5 and u = 0.5; Panel (d): the 
dashed, solid and dotted lines correspond to u = 0, 0.5 and 1 respectively with fixed m = 10, a. = 0.5 and a = 0.15. 
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= 
e 
o 


EN 
o 


o 
Rut 


Power-law flux/(10"?erg:cm ?-s^!) 


10; 


Power-law flux/(10 *erg:cm ?-s^!) 


= 
[2] 


E 


eo 
a 
eo 


0.01 


Power-law flux/(10 9 erg:cm ?:s^!) 


0.01 


Type-B QPO frequency/Hz 


Power-law flux/(10 9 erg-cm ?:s^!) 
a 


Type-B QPO frequency/Hz 


Fig. 2 The theoretical QPO-PL relation with i = 50° and d = 5 kpc. The fixed parameters are adopted: (a) a, = 0.5, a = 0.15, 


u = 0.5; (b) m 


10, a= 0.15, u = 0.5; (c) m 


10, a, 


0.5, p = 0.5; (d) m = 10, a, 0.15. 


In Fig. 3, we compare the modelled correla- 
tion with the observed one (Fig. 8 in Ref. [34]). 
We adjust the values of a and p to fit the observed 
relation with fixed m = 10, a, = 0.5, i = 50? and 
d — 5 kpc. The modelled QPO-PL relations fit well 
into the observed three subgroups with a ~ 0.1- 
0.3 and u ~ 0.2-0.4. A smaller a value and a larger 
p value correspond to the upper subgroup. 

As shown in Fig. 8 of Ref. [34], the observed 
QPO frequencies and the power-law flux for three 
single sources, namely, GX 339-4, H1743-322 and 
XTE J18594-226, are well correlated. We fitted the 
data of these three sources separately as shown in 
Fig. 4, adopting the source parameters (BH mass, 
inclination and distance) listed in Tab. 1 of Ref. 
[34]. For H1743-322, We assume m — 10 since its 
BH mass is unknown and a, = 0.2/9 is taken. For 
GX 339-4, we adopt i — 50? following Ref. [47] 
and a, = 0.9448. For XTE J1859+226, we take 
a, = 0.5 since its BH spin is unknown. About the 
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0.5, a 


same value ps ~ 0.5 is needed to fit well into the ob- 


served correlation of these three sources as shown 


in Fig. 4. 
50 m~ T T 
E == a@=0.07,u=0.35 we 
T — @=0.15,u=0.4 AA 
o == a-0.34,u-0.2 A 
5 
E | 
9 depts 
2 10 O7 0x | 
7 d 
S5 ] 
E1 
5 * H1743-322 
= n XTE J1859«226 
2 „5 o GX 339-4 
T ^ x 4U 1543-47 
aL E A XTE J1550-564 
2 " n * XTE J1752-223 | | 
$ pe v XTE J1817-330 
a 并 + MAXI J1659-152 
0.5 oe oq EE pr borob ok r lio 3 T 1 E X 1 1 I3 
3 4 5 6 7 8 


Type-B QPO frequency/Hz 


Fig. 3 Comparison of the theoretical QPO-PL relation with 
the observed QPO frequency and power-law flux for a 
sample of eight BHBs. The curves are plotted with m = 10, 
ax = 0.5, i = 50° and d= 5 kpc. The observational data 
(taken from Ref. [34], without error bars) of the eight BHBs 


are shown with different symbols. 
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Fig. 4 The QPO-PL relation of GX 339-4, H1743-322 and 
X'TE J18594-226. The curves represent the theoretical 
relation and different symbols represent the observational 


data as same as those in Fig. 3. 


4 Conclusion 


We explain the QPO-PL relation quantita- 
tively based on the MHD Alfven wave oscilla- 
tion model. The type-B QPO is produced by the 
Alfven wave oscillation in the inner accretion disk. 
The power-law flux is calculated based on inverse 
Compton scattering of soft photons from SSD by 
a medium of hot electrons. We find that the mod- 
elled correlation fits well into the observed one 


within reasonable parameter range. The results 
suggest that magnetic fields in the inner region 
of the accretion disk play a key role in under- 
standing type-B QPO and the associated power- 
law flux. The hot electrons could be heated by 
the magnetic reconnection of the tangled small- 
scale magnetic fields in the inner disk? 99], More 
energy is released in the magnetic reconnection 
with a stronger magnetic field and hence brings 
about a higher electron temperature and a stronger 
power-law flux. The positive correlation between 
the QPO frequency and power-law flux is expected. 

'The type-B QPO data used in this paper was 
observed in the rising phase of outbursts given in 
Ref. [34]. Data in the decay phase is not included. 
In order to test whether our model is applica- 
ble to type-B QPOs in the decay phase, we also 
compared our modelled QPO-PL relation with the 
observed one by Ref. [14] including data both in 
the rising and decay phase of GX 339-4. How- 
ever, the modelled relation does not agree with 
that observed in the decay phase. In the decay 
phase, both the observed power-law flux and QPO 
frequency are low. The slope of the correlation is 
gentle and does not become steep at low accretion 
rates as in our model (Figs. 2 and 3). In the de- 
cay phase, the accretion flow may be different, i.e., 
with different radiative efficiency, or the corona be- 
comes weak. The association of type-B QPOs and 
jet is supported by more and more observational 
evidence!?>*3], We noticed that the toroidal Alfven 
wave oscillation could also occur in the jet since the 
toroidal magnetic field will grow with the rotation 
of the BH (Blandford-Znajek! process) or the 
disk (Blandford-Payne!*?! process). The electrons 
in jet may be heated by the magnetic energy to 
produce the power-law flux. Therefore, our model 
is also applicable to the jet origin of LFQPOs. Fur- 


ther investigations are expected in the future work. 
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ale) Me HY BAY QPOs [Win ae B S BOTH TERJ 
ME 
TAN KY 


(长 江 大 学 物理 与 光电 工程 学 院 EUN 434023) 


摘要 观测 表明 ， 黑洞 双星 的 B 型 准 周期 振荡 (Quasi-Periodic Oscillation, QPO) 43€ 5; 5&8 Z [IR] ££ 4E IETH2 TE. 试图 
基于 阿尔 文 波 振荡 模型 定量 解释 该 相关 性 . 标准 薄 吸 积 盘 辐射 通 量 极 大 值 处 的 阿尔 文 波 振荡 产生 QPO. 标准 薄 盘 上 的 软 光 子 
与 晃 或 喷 流 基部 的 热电 子 介 质 发 生 逆 康 普 顿 散射 产生 守 律 通 量 . 通过 吸 积 率 的 连续 变化 , 得 到 QPO 频 率 与 窜 律 通 量 关 系 的 分 
析 解 和 数值 解 . 模拟 得 到 的 相关 性 在 合理 的 参数 范围 内 与 观测 值 相 吻合 . QPO 频 率 与 完 律 通 量 的 正 相关 性 可 以 理解 为 , 较 强 
的 磁场 导致 较 高 的 阿尔 文 波 频 率 和 较 高 的 电子 温度 从 而 得 到 较 高 的 容 律 通 量 . 结果 表明 B 型 QPO 可 能 与 吸 积 盘 或 喷 流 中 的 
环 向 磁场 的 活动 有 关 . 


关键 词 AR, MAR, 磁场 , 黑洞 物理 , 恒星 : 黑洞 , X 射 线 : 双星 
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